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ABSTRACT 


The  possibility  of  expressing  any  B-spline  as  a  positive  combination  of 
B-splines  on  a  finer  knot  sequence  is  used  to  give  a  simple  proof  of  the  total 
positivity  of  the  spline  collocation  matrix. 
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SIGNIFICANCE  AND  EXPLANATION 

\ 

The  total  positivity  of  the  spline  collocation  matrix  Is  the  basis  of 
several  important  results  in  univariate  spline  theory.  This  makes  it 

desirable  to  provide  as  simple  as  possible  a  proof  of  this  total  positivity. 

v :  s 

The  proofs  available  in  the  literature  don't  qualify  since  they- all  rely  on 

f  f  c  t  ^  r 

certain  determinant  identities  which  are  not  exactly  intuitive.  give  here 
a  proof  that  uses  nothing  more  than  Cramer's  rule  (hard  to  avoid  since  total 
positivity  is  a  statement  about  determinants)  and  the  geometrically  obvious 
fact  that  a  B-spline  can  always  be  written  as  a  positive  combination  of 
B-splines  on  a  finer  knot  sequence. 

Hie  geometric  intuition  appealed  to  here  stems  from  the  area  of  Computer- 
Aided  Design  in  which  a  spline  is  constructed  and  manipulated  through  its 
B-polygon,  a  broken  line  whose  vertices  correspond  to  the  B-spline 
coefficients.  If  a  knot  is  added  (to  provide  greater  potential  flexibility), 
the  new  B-polygon  is  obtained  by  interpolation  to  the  old.  This  has  led  Lane 
and  Riesenfeld  to  a  proof  of  the  variation  diminishing  property  of  the  spline 
collocation  matrix  and  is  shown  here  to  provide  a  proof  of  the  total 

positivity  as  well. 

/V 

\ 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MFC,  and  not  with  the  authors  of  this  report. 


A  GEOMETRIC  PROOF  OF  TOTAL  POSITIVITY  FOR  SPLINE  INTERPOLATION 
C.  de  Boor1  and  R.  DeVore1'2 

fl.  Introduction.  Perhaps  a  better  title  would  be  "Adding  a  knot  can  be 
illuminating"  since  the  purpose  of  this  note  is  to  show  how  this  idea  can  be  used  to  give 
simple  proofs  of  several  important  properties  of  B-splines,  including  the  total  positivity 
of  the  B-spline  collocation  matrix  and  the  sign  variation  diminishing  property  of  the  B- 
spline  representation.  He  show  that  variation  diminution  follows  immediately  from  the 
fact  that  a  B-spline  on  a  given  grid  is  a  non-negative  linear  combination  of  B-splines  on 
a  refined  grid.  He  use  the  same  fact  to  prove  the  non-negativity  of  any  minor  of  the 
collocation  matrix  and.  with  a  bit  more  care,  even  characterize  which  of  these  minors  are 
positive. 

The  total  positivity  of  the  collocation  matrix  was  originally  proved  by  S.  Karlin  [5] 
in  his  development  of  the  general  theory  of  total  positivity.  Later  C.  de  Boor  gave  a 
spline  specific  proof  (3].  In  both  cases,  variation  diminution  was  derived  as  a 
consequence  of  total  positivity.  He  obtain  both  properties  directly.  This  was  motivated 
in  part  by  the  work  of  J.  Lane  and  R.  Riesenfeld  (6] ,  who  gave  a  direct  proof  of  variation 
diminution  based  on  spline  evaluation  algorithms  used  in  computer-aided  design  which  can 
be  interpreted  as  "adding  knots".  But  we  follow  BShm's  idea  [1]  of  adding  one  knot  at  a 
time.  He  note  that  Jia  [4]  has  done  related  work  concerning  the  total  positivity  of  the 
discrete  B-spline  collocation  matrix. 

Let  k  >  0  be  a  fixed  integer  which  is  the  order  of  the  splines.  He  call 
t^  (ti)"+k  a  knot  sequence  if  t^  <  ti+1,  1  <  i  <  n  t  k  and  tj^  <  ti+k,  i  “  1,...,n. 
The  B-splines  of  order  k  for  this  knot  sequence  t^  are  given  by 
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(1.1)  Mt(x)  «-  Ml#t(x)  «-  <t1+k  -  ti)(ti....,ti^(](*  -  *i'-\  i  -  1,...,n  , 
vhsrs  dtnotii  a  k-th  order  divided  difference  and  i"  auuc{u,0).  It 

follows  that  >  0  and  aupp  M^  “  ^ti>ti4k^*  0,1  M°h  interval  (tj.tj Mj  ia  a 
polynomial  of  order  k  (degree  <  k).  The  B-aplinea  are  linearly  independent  and 
|  Hj  :  1  on  [tk,tn].  In  particular  if  the  n unbar  x  •  (t1(tn4k)  appeara  exactly 
k-1  tinea  in  t,  then  there  ia  only  one  B-apline  which  ia  non-sero  at  x  and  ita  value 
at  x  ia  one.  For  thaae  and  other  propertied  of  B-aplinea,  see  [3]. 


$2.  Knot  refinenent.  Me  nay  that  the  knot  sequence  *  is  a  ref  I  ewe  at  of  t  if 
a.  contains  t  as  a  subsequence.  Our  only  tool  in  the  subsequent  arguments  is  the 
observation  that 

(2.1)  any  B-apline  -  Mj>t  ia  a  positive  linear  coafolnation  of  sow  of  the 
B-aplinea  HJ  »■  for  the  refined  knot  sequence  a.  .  Precisely, 

"j  -  i  0j(i>  »• 

with  dj  nonneqatlve,  and  aupp  -  [t,u]  ,  where  (at'*u+k)  is  the  snsllest 
swt  of  b  containing  (tj, . .  .,t  j*.k)  as  a  subsequence. 

Me  first  prove  (2.1)  for  the  special  case  that 


£  *  (  •••#  ty.J,  Sy,  ty,  ...  )  , 

i.e.,  B_  is  obtained  fron  _t  by  the  addition  of  the  knot  ay  (satisfying  tv_1  <  sv  < 
tv  ,  of  course).  Then 


(2.2) 


J 


M*  for  j4k  <  v 
for  v  <  J 


For  j  <  v  <  j+k  ,  we  have  two  ways  of  writing  the  divided  difference  (s^. 


j' 


j  . 

***'  *j+k+1-y 


■j+k+1  ■) 

with  S4  j»  (s1,...,s1+k)  ,  Tt  i-  [ti,...,ti+k]  .  Therefore 

<fcj+k  "  tj,Ti  ”  <«v  "  “i*8)  +  ^jek+l  -  «v>8j*1 


hence 


-2- 


(2.3) 


N  , 


SJ  ♦ 


j-flt-M  .v  jj  ,  j  <  v  <  j+k  . 


"  j+k  "j  J  j+k+1  "J+1 

He  can  combine  this  with  (2.2)  into  one  formula,  as  follows! 


(2.4a) 

with 

(2.4b) 

Consequently , 


"j 


(1 


Yj)Nj  +  Yj+1Nj+1  '  a11  j  ' 


min  { 


(« 


J1L 


*u)+ 


Vk  -  *3 


,  1  }  ,  all  j  . 


tj) 

,  if 

1  j+-k  *  av 

(2.5) 

supp  Uj  ■  [j»j+1) 

,  if 

<  -v  <  fc: 

tj+1] 

,  if 

8  <  t 

V 

j 

and  this  finishes  the  proof  of  (2.1)  for  this  case. 

The  general  case  follows  from  the  repeated  application  of  this  special  case,  by 
induction!  Suppose  that  £  is,  in  turn,  a  refinement  of  s.  ,  hence 

»»■  -  l  NJ  , 

with  Nj  1-  •  Then  it  follows  that 

(2.6)  Nj  ”  ^  8j(t)  "Z  '  With  6j<4>  "  ^i  aj{i)oi,t)  ’ 

Therefore  6j  >  0  since  we  already  know  that  Oj  ,  >  0  .  Further 


supp  0j 


iSauppu. 


U  supp  o| 


[i',u')  , 


j 


with  (rti,  ...»  ru<+k>  the  smallest  segment  of  r  containing  (S(,  •••,  «,*«*)  ••  a 

subsequence.  But,  since  H.u]  is  the  support  of  Uj  ,  i.e.,  (Sj,  ...»  su+k)  is  the 
smallest  segsMnt  of  _e  containing  (tj,...,tj+k)  ,  it  follows  that  (rt> , . . .  ,ru.+){)  is 
also  the  ssiallest  segment  of  r_  containing  (tj, •  •  •  ,tj+y)  • 

The  coefficient  function  Oj  in  (2.1)  has  been  called  a  discrete  B- spline.  The 
above  argument  shows  that  the  matrix  (a^ti)]  ie  the  product  of  bi -diagonal  matrices  with 
nonnegative  entries,  hence  totally  positive  by  Cauchy-Binet.  This  is  the  basic  idea  behind 
the  proof  of  such  total  positivity  in  Jia  [4] . 
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■,  i 

i  .1 

:  J 

!  J 

i  ' 

i  -i 

i  i 

A 

•  * 


1 


|3.  Variation  diminution.  He  use  tha  customary  notation  8“( a)  for  tha  number  of 

(strong)  aign  changes  in  tha  sequence  or  function  a.  Ha  want  to  shew  that  S"(E  XjNj)  < 

n 

S“(JO,  i.e.,  tha  spline  f  i*  [  X.H  changaa  aign  no  more  oftan  than  ita  coefficient 

,33 

sequence  X.  This  follows  frost 

(3.1)  1)  If  t  tm  )  XjHj  -  l  X^HJ  with  N j  i-  and  £  a  refinement  of  £  ,  than 

S-(X)  >S*(XM  » 

11)  If.  In  addition,  x  *  (t^.t,,**)  appears  as  a  knot  In  £  with  (exact) 
aailtlplldty  k-1  ,  then  X^  «  f(x)  for  aoss  j  . 

Property  11)  is  clear.  To  prove  property  1),  we  first  consider  tha  special  case 
when  £  is  obtained  fros  t  by  tha  addition  of  a  single  knot.  In  that  case,  wa  infer 
fros  (2.4)  that 

i  XjHj  -  l  X,((1  -  VjM*}  ♦  TJ+1HJ+1)  . 

Therefore 

(3.2)  l  XjHj  -  l  XJHJ  with  X*  TjXj.,  ♦  d-Y^Xj  .  all  3  • 

(Hera,  we  sat  X„  i-  0  .)  Since  8  (0,1],  this  laplles  that  S'U^.Xj.Xj)  * 

8  OjwpXj)i  Therefore  8  (^1)  •  S  (  «..,X^W,,  X^,  X^,  Xj^,i  ...  )  )  8  (j)1 )  . 

This  shows  (3.1.1)  for  a  single  knot  refinesent.  But  then  by  induction  wa  get  (3.1.1)  for 
any  refinesent- 

Theorem  1.  (Variation  Diminishing  Property).  S"( E  XjHj)  *  8” (.X)  . 

Proof.  Let  f  ■  2,"  X^Mj  .  He  want  to  show  that,  for  any  increasing  real  sequence 
(sjf  ,  8“(  (f  (*,) ))  <  8™(_X)  .  He  can  assume  that  the  z^  are  not  knots  and  that  z,  6 
(ti,tR+k*  (since  f  2  0  outside  this  interval).  Let  £  be  a  knot  refineannt  of  £  such 
that  each  s,  appears  exactly  k-1  times  in  s.  Then  from  (3.1.11)  the  sequence 
(f(sl>)‘  is  a  subsequence  of  X  and  the  dealrad  result  follows  from  (3.1.1).  ||| 

It  is  aoamtlmea  useful  to  visualise  the  coefficients  (X^)  geometrically.  If 
tj  «»  (tj*,  ♦  •••  ♦  tj>jt_,)/(k  -  1),  then  the  continuous  piecewlae  linear  function 
P(f»t.)  with  vertices  (t^.X^),  j  »  1 


n  is  called  the  B-polygom  of  f.  This  polygon 


changes  sign  exactly  as  often  as  For  a  single  knot  refinement  s_  of  t,  the 

points  s^  are  related  to  t*  as  in  (3.2),  i.e., 

sj  =  *1*1-1  +  (1-Vfcj  • 

Hence  the  vertices  of  P(f,sj  lie  on  P(f,t^;  which  is  another  way  of  viewing  property 
(3.1.1). 


§4.  Spline  interpolation.  We  now  consider  spline  interpolation  at  nodes 
Xi  <  xj  <  •••  <  xn  (later  we  allow  coalescence).  Given  (yj >",  we  have  the 
interpolation  problem 


(4.1) 


1  =  yi( 

j-i  3  3 


i  -  1, . . .  ,n 


with  coefficient  matrix 

(4.2)  A  At  (N j (*i ) )", j_i  • 

In  case  xA  »  tj,  we  require  that  this  point  appear  at  most  a  total  of  k  times  in  x 
and  t. 

We  will  show  that  A  is  totally  positive  and  furthermore  characterize  which  minors 
of  A  are  strictly  positive.  For  this,  let  B  be  a  square  submatrix  of  A  , 

B  -  A(I,J>  t.NjCxi))ieI(jej  , 

with  I  and  J  subsequences  of  (1,2,..-,n)  of  the  same  length. 


I  (i,,  ....  im)  ,  J  =*:  <jv  / 

say.  We  call  such  a  submatrix  "good"  if  all  its  diagonal  entries  are  nonzero.  This  is  a 

natural  distinction  to  make  here  because 

(4.3)  _i£  B  la  not  "good",  then  det  B  *  0  . 


Indeed,  assume  that  N,  (x.  )  »  0  for  some  p  .  Then  x,  0  (t.  ,t.  .v)  •  Assume  that 
Jp  ip  ip  3p  Ip-"* 


x^  <  tj  .  Then  Nj(Xg)  -  0  for  q  <  ip  ,  j  >  jp  ,  and  this  shows  that  columns  p, 


. ..,  m  of  B  have  nonzero  entries  only  in  rows  p+1,  ...»  m  ,  hence  are  linearly 


_+k 


dependent.  So,  det  B  *  0 


The  argument  for  the  case  x. 


is  similar. 


Next,  we  write  det  B  as  a  linear  coablnetion  of  determinants  of  the  form  A'(I,K) 


with 

A’  «-  (N'(xi)) 

and  (HJ)  the  B-splines  for  a  refinement  a_  of  t  .  Precisely,  we  claim  that,  for  e 
certain  nonnegative  aj  , 

(4.4a)  det  A(I,J>  -  l* *j(K)  det  A*(I,K) 

with  the  euperecript  "+*  indicating  that  the  eum  ia  only  over  lmcraaalmg  K  .  Further. 

(4.4b)  supp  aj  _  aupp  Oj  , 

where  Oj(K)  :•  a^tk^)  •••  <»j  (kg)  and  the  are  aa  in  (2.1). 

For  the  proof,  we  coneider  firat  the  apecial  case  that  a,  ia  obtained  from  t_  by 

the  addition  of  a  aingle  knot.  Since  Nj  ■  \  o^(i)NJ  by  (2.1),  the  linearity  of  the 

determinant  aa  a  function  of  the  columne  given 

(4.5)  det  A( I,  J)  -  l  OjUC)  det  A*  (I,K) 

with  aj(K)  »■  a^(kj)  •••  Oj  (kg)  .  Bacall  from  (2.5)  that  aupp  <>j  C  [J,j+1]  . 

Therefore,  retaining  in  (4.5)  only  terms  with  a^CK)  f  0  ,  we  have  kp  ■  jp  or  jp+1  , 

all  p  .  Thus  K  ia  strictly  increasing  unless  kp  -  kp+1  for  acme  p  (possible  in 

case  jp+1  -  jp^f  ).  But  in  the  latter  case,  the  determinant  ia  trivially  zero  and  hence 

can  be  ignored.  Thin  finishes  the  proof  of  (4.4)  for  this  special  case. 

He  prove  the  general  case  by  induction  an  the  length  difference  d  »■  |sj  -  |tj  , 

having  just  proved  it  for  d  -  1  .  Assuming  it  correct  for  a  given  d  ,  let  t  be  a 

refinement  of  t,  with  |rj  -  |tj  •  d+1  and  let  i  be  t  one-point  refinement  of  |t] 

which  is  refined  by  .  Then,  with 

A"  I-  (n?(x. ))  and  NJ  x-  N .  ,  all  j  , 

3  1  3  3*£ 

we  have  NJ  •  Z  aj(£)  Nj  •  Further,  fro®  (4*5)  and  the  induction  hypothesis, 

det  A(I,J)  -  i+bj(L)  det  A*(I,L) 

with 

(4.6)  bj(L)  s-  I*  dj(K)a£(I,)  >  0  , 

which  makes  (4.4a)  obvious. 


The  proof  of  (4.4b)  is  a  bit  more  complicated.  It  can  be  skipped  if  only  the  total 


positivity  of  A  is  of  interest.  We  must  show  that  supp  bj  ■  supp  Bj  .  with  Bj(L)  :* 
•jJV  ”•  8j.'V  *  Suppose  first  that  Sj(L)  =  0  .  Then  Bjd)  =  0  for  some  j  6  J  , 
1  0  L  .  Therefore,  from  (2.6),  £  <jj(i)  a[(t)  =  0  ,  and,  since  all  terms  in  this  sum  are 
nonnegative,  they  must  all  be  zero.  Thus,  aj(X)c£(I.)  »  0  for  all  K  .  But  by  induction 
hypothesis,  supp  a£  -  supp  ,  therefore  also  <Xj(K)  a£(L)  *  0  for  all  X  .  We 
conclude  with  (4.6)  that  supp  bj  C  supp  Bj  • 

To  see  that  supp  bj  O  supp  Bj  ,  we  must  show  that 

(4.7)  Bj(L)  9  0  implies  Oj(X)  o£(L)  +  0  for  some  increasing  K  . 

Since  supp  a£  *  supp  ,  this  implies  that  a(K)a£(L)  /  0  for  this  increasing  K  , 
hence  also  bj(I.)  ft  0  from  (4.6). 

For  the  proof  of  (4.7),  it  is  sufficient  to  show  the  existence  of  a  K  with 

(4.8)  k  e  A  {i  :  o.  (i)o!(l  )  f  0>  ,  all  p  , 

v  Jp  1  p 

and  kp  <  kp+1  ,  all  p  .  Since 

BjU)  -  ajtjjajd)  +  flLj(j+1)aJ+1U)  , 

6j(I>)  ft  0  implies  that 

9  9  Aj  C  {j,j+1}  ,  all  j  6  J  . 

Renee,  the  existence  of  X  satisfying  (4.8)  is  assured.  To  finish  the  proof,  we  must 

show  that  it  is  possible  to  choose  such  a  X  which  is  also  increasing.  If  A,  f)A,  = 

3P  Jp+1 

9  ,  then  we  have  kp  <  kp+1  for  any  X  satisfying  (4.8).  Thus  we  only  have  to  consider 

how  to  choose  the  components  of  x  corresponding  to  a  connected  component  A,  ,  . . . , 

Jp 

A.  .  By  this  we  mean  that 
•*<1 

n  Ai  *  9  for  P  <  V  <  q  , 

Jv  Jv+1 

while,  for  any  i  +  jp,  ...,  jq  , 

A,  0  1,  "  9  for  p  4  V  <  q  . 

Jv 

Then  we  can  write  (jp,...,jq)  “  ( j, j+1, • • * , j ' )  ,  hence,  q-p  -  j'-j  .  Further,  i  e  Aj 
for  i  «  j+1  ,  ....  j'  .  Hence,  if  also  j  e  Aj  ,  then  the  choice  kv  «  jv  ,  all  v  ,  will 
do.  In  the  same  way,  we  have  i+1  e  AA  for  i  «  j,  ...,  j'-1  .  Hence,  if  j'+1  e  Aj.  , 
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then  the  choice  ky  -  jv+1  ,  all  w  ,  will  do.  We  claia  that  the  remaining  case 

j  9  Aj  and  j'+1  9  Ajt 

cannot  occur  since  it  would  imply  that  there  are  at  least  k  entries  in  t_  between 
rt  and  rtptk  •  Indeed,  with  supp  aj  «s  (*,ul,  it  would  follow  that  u  <  ip  ,  while 
also  iq  <  i'  ,  with  supp  aj.+1  «i  ti'.u']  .  Further,  let  sv  be  the  additional  knot 
in  s  .  Then,  by  (2.5),  At  D  A1+1  t  9  implies  supp  -  U,i+1}  ,  hence,  by  (2.5),  si  < 
sv  <  si+k  ,  i“J,...,j*  ,  therefore  <  *j+k  *  and  80  *'  <  u+lc  w*111®  «l»o  p+k  -  q-1 

«  j+k  -  (j*+1)  <  u+k  -  t*  .  This  would  iaply  that 

*p  <  Vi  <  •••  <  »q  <  *'  <  u+k  <  *p+lc  ' 
hence  k  *  ip+k  -  (^  >  1  +  (u+k-f)  +  1  +  q-p  >  1  +  (p+k-q_1)  +  1  +  q-p  *  k+1  • 

Theorem  2.  The  mstrlx  A  of  (4.2)  is  totally  posltl'  Moreover,  the  aubmatrlx 

B  of  A  formed  by  rows  i^,...,^  and  columns  as  »  positive  determinant  if 

and  only  If  It  la  “good1*,  l.a., 

*  e  supp  h.  ,  v  •  l,...,i 
v  3W 

Proof.  We  already  proved  that  det  B  «  0  unless  B  is  "good".  Bow,  to  prove  that  a 
"good"  B  has  a  positive  determinant,  we  choose  a  refinement  ji  of  Jt  so  fine  that 

(4.9)  for  each  i  «  I  ,  Nj(x^)  9  0  implies  that  Nj(Xp)  ■  0  for  all  p  ^  i  . 

Then  each  A'(I,K)  appearing  in  (4.4a)  has  at  most  one  nonsero  entry  in  each  column, 

hence  is  "good",  therefore  nonsero,  only  if  it  is  diagonal,  in  which  case  its  determinant 

is  obviously  positive.  To  finish  the  proof,  we  must  show  that  at  least  one  of  the 

aatrices  appearing  in  the  Sum  in  (4.4a)  with  a  positive  coefficient  is  "good".  Here  is 

one  such.  Choose  K  so  that  sb  is  the  first  point  in  s  to  the  left  of  x,  , 

P  P 

p**1,...,a  .  Since  N,  (x,  )  F  0  ,  this  implies  that  a,  (k_)  y  0  ,  all  p  .  ||| 

Jp  p  Jp  * 

Corollary.  (I.  Schoenberg  and  A.  Whitney  [7J).  The  interpolation  problem  (4.1)  has 
a  unique  solution  for  all  (y^)"  if  and  only  if  x^  e  supp  N^,  i  -  1,...,n. 
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We  can  also  allow  coalescence  of  the  interpolation  nodes.  If  (z^)  is  such  a 
nondecreasing  sequence  of  nodes,  then  we  can  think  of  it  as  the  limit  of  strictly 
increasing  sequences  (x^)  .  Correspondingly,  repetition  of  a  z^  corresponds  to  repeated 
or  osculatory  interpolation,  i.e.,  the  matching  of  higher  derivatives.  Precisely,  (4.1) 
becomes 


(4.10) 


n  u. 

I  M  V(z.  )  -  yif 

1  J  J 


1 , . . .  ,  n 


where  is  the  number  of  j  <  i  for  which  Zj  =  z^.  We  still  require  that  any  point 

appear  at  most  k  times  totally  in  z  and  t.  The  coefficient  matrix  of  (4.7)  is 


(4.11) 


A 


A 


t  ,z_ 


(D 


vv> 


n 

i.j*1  ‘ 


It  is  clear  that  A  need  not  be  totally  positive  since  entries  involving  derivatives 
may  be  negative.  However,  as  a  well-known  argument  shows,  if  M  is  a  minor  formed  by 
rows  i1,...,im  and  columns  with  the  property 

(4.12)  i  ,  <  i  -  1  implies  z  <  zi  ,  v  =  1,...,m  , 

v"  1  v”  ’  1  v 

then  M  >  0.  In  fact,  if  M(jc)  denotes  a  minor  corresponding  to  distinct  nodes 
x^  ”  (xi,...,xm),  then  subtracting  row  one  from  row  two  shows  that  M(>t)/(X2  -  x^  ) 
converges  as  x2  *  x1  to  the  minor  M'  which  replaces  row  two  of  MOO  by  first 
derivatives  at  x^  Hence  M'  >  0.  Using  this  type  of  limiting  process  we  see  that  any 
minor  M  satisfying  (4.12)  is  non-negative. 

We  can  also  characterize  those  M  satisfying  (4.12)  which  are  positive,  namely,  they 
satisfy 

(4.13)  z^  e  supp  N.  ,  v  =  1,...,m  . 

v 

The  necessity  of  (4.13)  is  proved  in  the  same  way  that  the  necessity  of  (4.9)  was 
established. 


The  sufficiency  of  (4.13)  is  proved  by  making  slight  modifications  to  the  earlier 
proof.  For  this,  it  will  be  convenient  to  allow  a  point  z^  to  appear  a  total  of  more 


than  k  times  in  a  and  x.  This  ia  acceptable  provided  we  stipulate  that  all  B- 

eplinea  and  their  derivatives  be  interpreted  as  right  limits  at  such  ij,  that  is  at 

*  J.  with  this,  let  s.  be  a  refinement  of  £  such  that  each  node  s^  appears  as  a  knot 

in  s.  exactly  k  times,  and  similarly  sach  t^  appears  in  s.  exactly  k  times.  If 

J  satisfies  (4,13),  we  choose  L  so  that  a^  *  and  the  number  of  j  <  1  with 

P  P  P 

s,  -  s.  is  u,  •  Since  the  coefficients  o(K)  in  (4,4a)  are  independent  of  £  ,  we 
3  P  P 

then  obtain  det  A(I,J)  as  a  positive  combination  of  certain  (nonnegative)  minors  of  A* 

! „  Am  g  •  In  particular,  the  submatrix  A'(J,L)  will  appear  in  that  sum  with  positive 
coefficient  since  Oj(L)  >  0  ,  and  det  a'(J,l)  >  0  since  A'(J,L)  is  lower  triangular 
with  positive  diagonal.  He  have  therefore  proved  the  following  theorem. 


Theorem  3,  For  the  matrix  A  of.  (4.11),  and  each  I,  J  satisfying  (4,12), 
det  A(1,J)  >  0  .  This  minor  is  positive  if  and  only  if  (4.13)  is  satisfied.  In 
particular  (4.10)  has  a  unique  solution  if  and  only  if  e  supp  H^,  1  »  1,...,n. 
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